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Abstract 



Properties of the proton and the neutron are explored in partially- 
quenched chiral perturbation theory with two non-degenerate light flavors. 
Masses, magnetic moments, matrix elements of isovector twist-2 operators 
and axial-vector currents are computed at the one-loop level in the chiral 
expansion. 
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I. INTRODUCTION 



At some point in time lattice calculations will be able to determine the properties and 
interactions of the low- lying hadrons directly from QCD. Impressive progress is being made 
toward this goal, however present computational limitations necessitate the use of quark 
masses that are significantly larger than those of nature. Typically one has m|j!^**' ~ 500 MeV. 
In order to make a connection between lattice calculations of the foreseeable future and 
nature, an extrapolation in the quark masses is required. Chiral perturbation theory (xPT) 
provides a systematic description of low-energy QCD near the chiral limit, and this technique 
has been extended to describe both quenched QCD (QQCD) |T|-^ and partially-quenched 
QCD (PQQCD) The hope is that future lattice simulations can be performed with 

quark masses that are small enough to guarantee a convergent chiral expansion, thus allowing 
a meaningful extrapolation to the quark masses of nature. 

Recently, it has been shown how to include the low-lying octet and decuplet of baryons |^ 
into partially-quenched chiral perturbation theory (PQxPT) As the pion mass 

presently obtainable in lattice simulations is not much smaller than the physical kaon mass, 
it is natural to consider a theory with three light valence- quarks and three light sea-quarks. 
The masses, magnetic moments and matrix elements of isovector twist-2 operators were 
computed at the one-loop level in SU{6\3)l ® SU{6\3)r PQxPT [|]. The strange quark 
mass was taken to be different from the other two degenerate light quarks (isospin limit) 
and one of the sea-quarks was taken to have a mass different from the other two. However, 
many lattice simulations are performed with just two light quarks. In this work we deter- 
mine several observables of the proton and the neutron at one-loop level in PQxPT with 
two non-degenerate light quarks using SU{4:\2)l ® S'f/(4|2)/j PQxPT. Isospin violation of 
electromagnetic origin is not considered. 

In Section |I| we develop the general formalism for including baryons in SU{4:\2)l ® 
SU{4:\2)ii PQxPT. The masses and magnetic moments of the proton and neutron are com- 
puted to one-loop order in the chiral expansion in Section |ITT| and respectively. In 
Section and ^ we present one-loop level expressions for matrix elements of isovector 
twist-2 operators and axial- vector currents, respectively. We conclude in Section [VTI. 



II. PQxPT 

The Lagrange density describing the quark-sector of PQQCD is 

a,b=u,d a,b=u,d a,b=j,l 

= Q' [^JP-^q]1 Qn , (1) 

k,n=u,d,u,d,j,l 

where the qv are the two light-quarks, u and d, the q are two light bosonic (ghost) quarks 
u and (i, and the gsea are the two sea-quarks j and /. The left- and right-handed valence-, 
sea-, and ghost-quarks are combined into column vectors 

Ql = (u,dJJ,u,dj^ , Qr = {u,d,j,l,u,dj^ , (2) 
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where the graded equal-time commutation relations for two fields is 

Q?(x)Qf (y) - {-r''Qi\y)Qm = 6"'6,,5%^ - y) , (3) 

where are spin-indices and i,k are fiavor indices. The objects r^^ correspond to the 
parity of the component of Qk, with 77^, = +1 for k = 1,2,3,4 and rjk = for k = 5,6, 
and the graded commutation relations for two Q's or two Q^'s are analogous. The Ql,r in 
eq. transform in the fundamental representation of SU{4:\2)l^ji respectively. The ground 
floor of Ql transforms as a (4, 1) of SU{4:)qL ® SU{2)qL while the first fioor transforms as 
(1, 2), and the right handed field Qr transforms analogously. 

In the absence of quark masses, mq = 0, the Lagrange density in eq. ([|) has a graded 
symmetry U{4\2)l ® f/(4|2)/j, where the left- and right-handed quark fields transform as 
Ql ^ UlQl and Qr UrQr respectively. The strong anomaly reduces the symmetry of 
the theory, which can be taken to be SU (4|2)i(g)5'f/ {^2) r®U {l)v [0]- It is assumed that this 
symmetry is spontaneously broken according to the pattern SU {A\2)l®SU {A\2) r®U {1)v — >■ 
SU{A\2)v ® U{l)v so that an identification with QCD can be made. 

The mass-matrix, mq, has entries mq = diag(m„, m^, rrij, mi, m„, m^), i.e. m^ = m^ and 
mj = m^, so that the contribution to the determinant in the path integral from integrating 
over the qy and the q exactly cancel, leaving the contribution from the gsea alone. 



A. The Pseudo-Goldstone Bosons 



The strong interaction dynamics of the pseudo-Goldstone bosons are described at leading 
order (LO) in PQxPT by a Lagrange density of the form 0, 



C = — str 



S^St^^S ] + A ^str [ mgSt + m^S ] + a^d^^d^^o - ml^l , (4) 



where 0$ and mo are quantities that do not vanish in the chiral limit. The meson field is 
incorporated in S via 



exp 



2t<!>' 
. f , 



(M 



(5) 



where M and M are matrices containing bosonic mesons while x and are matrices con- 
taining fermionic mesons, with 



M 



( nu 


7r+ 


JO 






Vd 


J- 






J+ 






U~ 






Vi 1 



M 



X 



Xrju 
Xtt" 



X7r+ 

Xvd 



Xjo Xl+ 
XJ' Xl" 



, (6) 



where the upper 2x2 block of M is the usual triplet plus singlet of pseudo-scalar mesons while 
the remaining entries correspond to mesons formed with the sea-quarks. The convention we 
use corresponds to / ~ 132 MeV, and the charge assignments have been made using an 
electromagnetic charge matrix, Q^^'^) = |diag (2, —1, 2, —1, 2, —1). 
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The singlet field is defined to be $o = str ( $ ) / y^, and its mass mo can be taken to 



be of order the scale of chiral symmetry breaking, niQ A. 



In taking this limit, one 



finds that the 77 two-point functions deviate from simple single-poles. The rjarib propagator 
is found to be 



id" 



{q^ - m]^){q^ - ml 



mt 



it 2 (g2 — m^^ + it){q'^ — m^^ + it){q'^ 



+ ie) 



(7) 



where mab is the mass of the meson composed of (anti)-quarks of flavor a and h. This can 
be compactly written as 



S^'Pa + nab{Pa,Pb,Px) 



(8) 



where 



P 



Hab{A,B,C) 



2 ^2 \(^2 ^2 



P. 



X 



q^ — rrix 



q^ — m^^ + le 
1 



A 



B 



+ 



{m\ - m]^){m\ - ml) 
[m\ - m^^J{m\ - m^^) 



C 



(9) 



where the mass, m^, is given by m\ = ^ (jT^jj + One important implication of this 

form for the singlet propagator is that the nucleon-nucleon potential falls off exponentially 
at long distances, as opposed to Yukawa-like behavior, away from the QCD-limit |^. 
The tree-level masses resulting from the Lagrange density in eq. (H) are 



m. 



ud 



m 



A {mu + md) , m. 



2 A mu 



(10) 



and so forth. 



B. The Nucleons and A- Resonances 



It is straightforward to include the proton, neutron, and A- resonances into PQxPT . 
An interpolating field that has non-zero overlap with the nucleons (when the ijk indices are 
restricted to 1,2) is [ffl 



ijk 



(11) 



where C is the charge conjugation operator, a, 6, c are color indices and a,/5,7 are Dirac 
indices. Dropping the Dirac index, one finds that under the interchange of flavor indices [ffl , 







(12) 



In analogy with QCD, we consider the transformation of Bijk under SU{A\2)y transforma- 
tions, and using the graded relation 
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Q^ u\ = u\ , (13) 

in eq. (|l3); it is straightforward to show that 

where Bijk describes a 70 dimensional representation of S'[/(4|2)y. 

It is convenient to decompose the irreducible representations of SU{4\2)v into irreducible 
representations of SU (2)vai (S> SU (2)sea ® SU (2)^ p,|TO|-Q. The subscript denotes where the 
SU{2) acts, either on the qv, the gsea, or the q. In order to locate a particular baryon 
in the irreducible representation we employ the terminology of Ref. ground floor, first 
floor, second floor and so on, as it is common in the description of super-algebra multiplets. 
The ground floor contains all the baryons that do not contain a bosonic quark, the first 
floor contains all baryons that contain one bosonic quark, the second floor contains all 
baryons that contain two bosonic quarks, and the third floor contains the baryons that are 
comprised entirely of bosonic quarks. As a way of distinguishing between baryons containing 
some number of valence and sea quarks, we use "levels" [§]. Level A consists of baryons 
containing no sea quarks, level B consists of baryons containing one sea quarks, level C 
consists of baryons containing two sea quarks, and level D consists of baryons composed 
only of sea quarks. 

The ground floor of level A of the 70- dimensional representation contains nucleons, 
comprised of three valence quarks, qyqvlv, and is therefore a (2, 1, 1) of SU (2)vai®5't/(2)sea® 
SU{2)q. The nucleons are embedded as 

Babe = ^{eabN, + eacNk ) , (15) 

where the indices are restricted to take the values a,b,c = 1,2 only. The nucleon doublet is 

N-{1) . (16) 

The first floor of level A of the 70-dimensional representation contains baryons that are 
composed of two valence quarks and one ghost-quark, qqvqv, and therefore transforms as 
(3, 1, 2) © (1, 1, 2) of S'f/(2)vai ® 5'[/(2)sea ® SU{2)q . The tensor representation aSab of the 
(3, 1, 2) multiplet, and gt of the (1, 1, 2) multiplet, where a = 1, 2 runs over the q indices 
and a, 6 = 1, 2 run over the qy indices, have baryon assignments 

The right superscript denotes the third component of g^-isospin, while the left subscript 
denotes the q flavor. The ground floor of level B of the 70-dimensional representation 
contains baryons that are composed of two valence quarks and one sea quark, qvQvQsen, and 
therefore transforms as (3, 2, 1) © (1, 2, 1) of SU{2)^^i 5'[/(2)sea ® 5'[/(2)^ . The tensor 
representation aSbc of the (3, 2, 1) multiplet, and at of the (1, 2, 1) multiplet, where a = 1, 2 
runs over the gsea indices and b,c = 1,2 run over the qv indices, have baryon assignment 

aS 11 = , aSl2 = aS21 = "^^a ' '''^22 = ^ at = A.^ ■ (18) 
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The aSab, at, aSab, 8Jid at 8.16 Uniquely embedded into Bijk (up to field redefinitions), 
constrained by the relations in eq. (|T^: 

Bijk = ~Y ^ i-2Sjk for I = 3,4 and j,k = 1,2 

1 1 

= 2 j-2t £ik + j~2Sik for j = 3,4 and i,k= 1,2 

Bijk = ^ k-2t Sij + for A; = 3, 4 and ii,j = l,2 

-Bijfc = W — i-4Sjfc for « = 5, 6 and j, A; = 1, 2 

%fc = ^ j~4^ £^ifc + j-4.Sik for j = 5, 6 and i,k= 1,2 

Bijk = -^ k-J Sij - k-4Sij for fc = 5,6 and i,j = l,2 . (19) 

We do not construct the remaining floors and levels of the 70 as we are only interested in 
one-loop contributions to observables with nucleons in the asymptotic states. 

As the mass splitting between the A-resonances and nucleons in QCD is much less than 
the scale of chiral symmetry breaking (A^ ~ 1 GeV) the A-resonances must be included as a 
dynamical field in order to have a theory where the natural scale of higher order interactions 
is set by A^^,. An interpolating field that contains the spin-| A-resonances as the ground 
floor is [ffl 



ijk ~ 



QrQf'Qr+Q^^'QrQr+QrQrQ^^' eabciCY)^, , (20) 



where the indices k run from 1 to 6. Neglecting Dirac indices, one finds that under the 
interchange of flavor indices 

= (-)'+''"*^-.fc = (-)^+'*'"=7;fc, . (21) 

Tijk describes a 44 dimensional representation of SU{A\2)y, which has the ground floor 
of level A transforming as (4, 1, 1) under S't/(2)vai ® S'f/(2)sea ® SU{2)q with Tabc = Tabc, 
when the indices are restricted to take the values a,b,c = 1,2, and where Tabc is the totally 
symmetric tensor containing the A-resonances, 

Tin = A++ , Tn2 = ^A+ , Tu2 = ^A° , T222 = A" . (22) 

The first floor of level A of the 44 transforms as a (3, 1, 2) under SU {2)^^i®SU{2)^^^(g>SU{2)g 
which has a tensor representation, aXij, with baryon assignment 

aXu = S^'^^ , aXu = aX2l = ' ^"^22 = ^ ■ (23) 

Similarly, the ground floor of level B of the 44 transforms as a (3, 2, 1) under S'f/(2)vai ® 
SU{2)sca ® SU{2)q which has a tensor representation, aXij, with baryon assignment 
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aXii — S*'^^ , aXl2 — aX2l — "^^a^ ' ""^22 — S*' ^ . (24) 

The embedding of aXij into Tijk is unique (up to field redefinitions), constrained by the 
symmetry properties in eq. (pl|): 



Tijk = +--^ i-2Xjk for i = 3,4: and j, k = 1,2 



1 
1 

73 



^ijk = -\ — 7^ i-2a;ifc for j = 3,4 and i,k = 1,2 



%jk = ^ — 7^ k~2Xij for fc = 3,4 and i,j = 1,2 



Tijk = H — 1= i-iXjk for i = 5,6 and j, k = 1,2 
v3 



T^jfc = --^ j-4a;iA; for i = 5,6 and /c = l,2 



T^jfc = H — 1= k-AXij for fc = 5,6 and i,j = 1,2 . (25) 
v3 

Again, we do not exphcitly construct the second and third fioor baryons of the 44. 



C. Lagrange Density for the Nucleons and A-Resonances 

The free Lagrange density for the Bijk and T^k fields is [Q, at LO in the heavy baryon 
expansion [p!^-[T7|, 



C = I (Bv ■ VB) + 2aZ^'^ (BBM+) + 2(3^j^^ (BM+B) + 2(t£*^^ (BB) str {M 



T'v . VT,) + A (tX) - 2jT,^^ {t'M^T,) - 2a'P {7%) str (M^) , (26) 



where A is the mass splitting between the 70 and the 44, Ai+ = (^'"^q^' + ^f^Q^j, and 

^ = y/r,. The brackets, ( ) denote contraction of Lorentz and fiavor indices as defined in 
Ref. [^]. For a matrix acting in spin-space, and a matrix Yij that acts in flavor-space, 
the required contractions are 

BTB) = B'^^'"' F^ , (t' F T„) = T'-^'^' 



{BTYB) = B"'"'" F^ Y, ' Bijk,p , ((T'TY T^) = T'^"''^' F^ Y, ' %^k,p, 

{BT by) = (-)(v^+V,){Vk+Vu)^'^^''^i y^n g..^^^ 

(b F F'^T^) = r F^ {Yn\ Ti,k,p, , (27) 
where B and T transform the same way, e.g. 

-gfcji _^ ^^■^^^(r}j+r}^)+{rii+ri^){rn,+rj„) -^nml ^ k\ jj j] jj^i] ^2g) 

The covariant derivative acting on either the B ot T fields has the form 
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where the vector and axial- vector meson fields are 



y' = \[ i^'i^ + i^Q'i ) ^ = \{ ^^^^^ - ^^^^^ ) 



(30) 



By restricting ourselves to the qvQvQv sector, we can make a comparison between the LO 
partially-quenched free Lagrange density and that of QCD, which has fewer free parameters, 



- iT^v -VT^ + A T% - 27mT'^A<2'^°T^ - 2aMT% tr 



M 



QCD 



(31) 



where the QCD chirally-invariant mass operator is Ai'^'^^ = | (■C^'^g^^ + 'C^q'C) j where 
rUg = diag(m„, m(i), and ^ is the QCD version of the matrix defined in eq. (H). The QCD 
parameters and the PQQCD parameters are related by. 



(PQ) 



{PQ) - -(PQ) 



(32) 



The Lagrange density describing the interactions of the 70 and 44 with the pseudo- 
Goldstone bosons at LO in the chiral expansion is [ffl. 



{Ta^b 



BA,r'' 



(33) 



where is the covariant spin- vector [|l3|-p!5[]. Restricting ourselves to the qvQvQv sector, 
we can compare eq. (Bl) with the LO interaction Lagrange density of QCD, 



C = 2gA NS'^A^N + g^NS^N tr [ A J + g^^ 
+ 2(7AA T'^^M^T, + 2gx TS^ tr [ A J , 



T ' A'^N.e,. + h.c 



and find that at tree-level 

4 1 
" = 7:9a + -9i , 



P 



^91 - ^9A 



n 



9aa 



C 



-9an 



(34) 



(35) 



with gx = 0. Considering only the nucleons, and decomposing the Lagrange density in 
eq. (|34D into the mass eigenstates of the isospin-symmetric limit, 7r^,7r° and rj, we have 



— ~ ^/2 — 

C = 2gANS^A^N + ^ {gA + gi) NS'^N d^r] , 



(36) 



where is the axial- vector field of pions only (excluding the isosinglet meson). In the 
isospin symmetric limit, with the mass of the 1] being of order ~ A^, all expressions must 
be independent of the coupling gi. 
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III. NUCLEON MASSES 



The mass of the i-th nucleon has a chiral expansion 

M, = Mo(/i) - M«(/i) - M^'/'\fi) + ... , (37) 

where a term M-"^ denotes a contribution of order ttlq, and i = p,n. The nucleon mass is 
dominated by a term in the PQxPT Lagrange density, Mq, that is independent of mg. Each 
of the contributions in eq. (0) depends upon the scale chosen to renormalize the theory. 
While at LO and next-to-leading order (NLO) in the chiral expansion, the objects Mq and 
M-^^ are scale independent, at one-loop level they are required to be scale dependent. The 
leading dependence upon uiq, occurring at order O^mg), is due to the terms in eq. (pHj) 
with coefficients o^m'^^ /^ir^^ and 

'^m'^\ 6ach of which must be determined from lattice sim- 
ulations. The leading non-analytic dependence upon mq arises from the one-loop diagrams 
shown in Fig. |I|, and we find the contributions to the proton mass are 




FIG. 1. One-loop graphs that give contributions of the form ~ mg to the masses of the proton 
and neutron. A solid, thick-solid and dashed line denote a 70-nucleon, 44-resonance, and a meson, 
respectively. The solid-squares denote an axial coupling. 



M(3/2) 
p 



9\ 



12 



+ 



9a9i 



m% + mfu - ^Id + 2ml^ + SCr,.^,^^ + 3G^„,^„ 



9an 

' 97T 



5Fud + Fuu + Fju + Flu + 2Fjd + 2Fici + 2i?^^_^^ + 2i?^^^^^ — AEr^^^r)^ 



(3J 



where the function F^, = F{mc, A, /x) is 



F{m, A, /i) = (m^ - A^) \^VA^ - log 
1 



'A - VA2 -m2 + ie^ 
A + VA2 - m2 + 



A log 



- -Am log ( - 



(39) 
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The functions G^^,^, and E"^^,^, are G^^^^j, 



'Hvavti^va^^vb^^x) and 



'^vaVbi^vay ^Vb^ ^x), respectively, where the function Hrj^rit is given in eq. (|). The con- 
tributions to the neutron mass are 



(PQ) ( I A 

\i {rrij + mi) 



9a 
3 



+ 



12 



"^dd - ^Kd + 2m-^d + ^^ju + 2mf^ + 3mf„ + 3G^,„,,„ + 6G,„,,, + 3G,,,^, 



9an 

' 9lT 



[ 5Fud + Fdd + Fjd + Fid + 2Fju + 2Fiu + 2Er,^^r]a + 2-E'r?„,r;u " 4-E,,„,^ 



(40) 



Our expressions for both the proton and neutron masses collapse down to those of isospin- 
symmetric QCD [|l^ in the limit mj,mi,mu,md — *■ m, 



Mn = Mo -2m{aM + 2cr 



M) 



^^2^3 I ^9 AN TT 



(41) 



In obtaining this result we have used the fact that Gn^^ria ~^ ~\^n^ and Erj^^riu ~^ By 
varying both the sea-quark and valence-quark masses over a suitable range and determining 
the proton and neutron masses, the constants Mq, o^m'^^ Pm'^^ and cr^'^^ can, in principle, 
be determined for two-flavor QCD. 



IV. NUCLEON MAGNETIC MOMENTS 

The nucleon magnetic moments have been computed on the lattice in both quenched 
and unquenched QCD using quark masses significantly larger than those of nature [0. In 
order to define the magnetic moments of the nucleons in PQQCD, one must first define the 
light-quark electric-charge matrix. As discussed in detail in Ref. [Q, the extension of the 
electric-charge matrix from QCD, where 

Q = diag(+^, , (42) 

to PQQCD is not unique, but is constrained by the requirement of recovering QCD in the 
limit that rrij rriu and mi md- The most general charge matrix whose matrix elements 
reduce down to those of QCD (keeping the valence-quark charges fixed) is 

= diag [+1,-1, 1j , qi, Qj , qi) ■ (43) 

For subsequent discussions we define Q[^^^ = ^ ( ^^Q^^'^^.^ + 'CQ^^'^^^^) i and Qj+ = 
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Studying the behavior of the magnetic moments for a variety of charge-matrices with 
different relative contributions of the singlet and adjoint representations will yield the desired 
low-energy constants in the chiral Lagrangian [p|,pO|. The additional freedom introduced by 
the charges qi and qj is a blessing rather than a curse due to the fact that different extensions 
correspond to different weightings of disconnected diagrams in lattice simulations. Thus the 
impact of disconnected diagrams whose numerical evaluation, in some cases, converges slowly 
(and hence can induce large uncertainties) can be minimized by a suitable choice of Q'^^^\ 
Furthermore, in PQxPT the contribution of one-loop diagrams involving the more massive 
sea-quarks can be minimized (order-by-order) by an appropriate choice of charges, thereby 
improving the convergence of the chiral expansion. 

In QCD, there are two invariants that can be constructed at LO in the chiral expansion 
to describe the magnetic moments of the nucleons, 

^ = Jir^^'' ( /^A Tr [ Q^^ ] Na^^N + /.^ Na^'QuN ) , (44) 



N 



where F^y is the electromagnetic field-strength tensor, and M^v is the physical value of the 
nucleon mass. More conventionally, this is written in terms of the isoscalar and isovector 
operators 

^ = JW^^'' ( ^"""^^ + iVa^V|+ N ) , (45) 



where /xq = (2/iA + fi'B)/G is the isoscalar nucleon magnetic moment, /ii = /iB/2 is the 
isovector nucleon magnetic moment, and r|Y = | ( ^^^'^^ + ^''""^^)- PQQCD, the LO 
Lagrange density contributing to the magnetic moments of the nucleons has the form 



am/"' 



+ /i^ str [ Qf/ 1 ( B a^'" B)] . (46) 



It is interesting to note that there is one more operator in PQQCD than there is in QCD at 
LO. The coefficients in the QCD and PQQCD LO Lagrange densities are related by 

/io = 7 [ + + ] , /il = i [ 2/i„ - /i^ ] . (47) 
o o 

At the order to which we will be working, my^, it is convenient to write the proton and 
neutron magnetic moments in terms of the /xo,i- However, higher order calculations will need 
to be performed using /ia,/?,^ directly. 

Up to order m^^, it is convenient to write the magnetic moment of the i-th nucleon as 

/^. = «^ + ^ [A + A'] + - , (48) 



where the tree-level contributions, a^, are given by the Lagrange density in eq. (46). The /? 
and (3' contributions arise from the one-loop diagrams shown in Fig. ^ involving nucleons in 
the 70 and resonances in the 44 dimensional representations, respectively. For the proton 
we find 
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N 



N 



N 



N 



T 



/ 

^ M 



T 



/ 

^ M 



FIG. 2. One-loop graphs that give contributions of the form ~ y/mg to t/ie magnetic moments 
of the nucleons. A solid, thick-solid and dashed line denote a 70-nucleon, 44-resonance, and a 



meson, respectively. The solid-squares denote an axial coupling from eq. (3S), while the solid circle 
denotes a minimal coupling to the electromagnetic field. 



ttp = /io + 1^1 

f3p = - - Aruju - 4:miu ] + [ mud + muu - rriju - rriiu 

ttt: [ T^ud + 4m„„ - Smrfrf + 3m,d - 4m,-„ + ?,mid - Ami 
36 



, /2^i r 1 , 9igA r 1 , ^1 " 



3 



+9« ( ^ [miu - mud] + [miu - + ^ 



3 3 
miu - ruud + -mid - 2"^<i<i 



/^p ~ 9\n [ -^dd ~ -^MU 6^ ud ~ ^ jd ~\~ ^ ju ~ 



6 

7d + -^in 



3 3 
+- g, ( + 2:Fud - - 2r,d) + - (-^.d + 2:Fdd - - ITi 

where the function Ti = J-'{mi, A, /i) 



T( A ^ /A2 21 f ^-V^'-m^ 

7ry-(m, A, u) = v A^ — m^log . = 



^ ~ ~ " ) - Alog(^ 

,A + VA2-m2 + ze; V / 

In the limit A ^ 0, JF(m, 0, /i) = m. For the neutron we find 



m 



2 

/3n = Y [7m„d + 2mid + 2mjd - 2mdd] + [mja + mid - mud - m^d] 

+1^ [3m„„ + 2mud - mad + m^d - ^mju + niM - 3m,„] 
io 



+<lj [mjd - mud] + [mjd - mud] + ^ 



mjd - mud + - o 



,2 r 



\ 3 3 

f^'n — 9aN 7T [ -^dd — 4JF + 9JF ud — ^ jd + 4JF — J-'id + 4:J^l 



3 3 
mid - mdd + -miu- -niud 



+3 Qj (JF + 2JF — ^ jd — 2JF + 3 qi {Tdd + 2JF — JF^^ — 



(49) 



(50) 



• (51) 
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These expressions reduce to their QCD counterparts pT|-p3 



fJ'p — f^o + /^i 



M, 



N 



47r/2 



2 2 'T~ 

9a "^71+ + g 9 AN -F-K^ 
2 



(52) 



when rrij — > and — m^. Note that the results are independent of the charges q^^i. 



V. MATRIX ELEMENTS OF ISOVECTOR TWIST-2 OPERATORS 

The forward matrix elements of twist-2 operators play an important role in hadronic 
structure as they are directly related to the moments of the parton distribution functions. 



There exist both quenched and unquenched lattice simulations of these matrix elements p4 
Recently, it was realized that the long-distance contributions to these matrix elements could 
be computed order-by-order in the chiral expansion using chiral perturbation theory ||25|-27 . 



These corrections have been applied to results from both quenched and unquenched lattice 
data p8|, with interesting results. In addition, the long-distance contributions arising in 



QQCD and the large-A'^c limit of QCD have been computed in Ref. p3 and Ref. [BD 



respectively. Further, this technique has been applied to the off-forward matrix elements 
of twist-2 operators in order to study the spin structure of the proton [^]. In Ref. |^, the 
matrix elements of the isovector twist-2 operators were computed at the one-loop level in 
S\J{^\?))l ® >S'f/(6|3)ii; PQxPT. Like the magnetic moments, the extension of the twist-2 
matrix elements from QCD to PQQCD introduces a non-unique isovector charge matrix, 
which can be exploited to optimize the numerical simulations and the chiral expansion 
In QCD, the nonsinglet operators have the form, 

^t},l ...M. = ^ ^{^^ - ^M.}) Q - traces , (53) 

where the {...} denotes symmetrization with respect to all Lorentz indices, and r" are Pauli 
matrices acting in flavor-space. They transform as (3, 1) © (1, 3) under SU{2)l ® SU{2)ji 
chiral transformations p5| , |2B| . Of particular interest to us are the isovector operators where 
= diag(l, —1). At LO in the chiral expansion the C^^j^jij ^„ match onto [p5| -p7| 




— traces . (54) 
In PQQCD, the nonsinglet operators have the form 

''^(^Si;: = ^Q^ B,,) ...{^ B,„}) Q - traces , (55) 
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where the r° are super Pauh matrices, an extension of the Pauh matrices from two-flavor 
QCD to PQQCD. With the requirements that is supertraceless and that the QCD matrix 
elements are recovered in the limit rrij — > mi — >■ m^, the most general flavor structure 
for is (keeping the valence-quark charges fixed) 



-1 , Vj , yi , Vj , yi) ■ (56) 



For an arbitrary choice of the yi, this operator contains both isovector and isoscalar com- 
ponents. It is purely isovector only when yj + yi = 0. As a result, for arbitrary y^, the 
usual isovector relations between matrix elements do not hold. The fact that disconnected 
diagrams can only be isoscalar renders this result obvious. At LO in the chiral expansion, 
matrix elements of the isovector operator ^'^C^^jjj ...^^ are reproduced by operators of the 
form ||,|5|-|3] 







t( 
















— traces 





1 \ 



str 



(57) 

where r|_,_ = | (■C'^'^'C^ + 'C^'^'^'C) • general, the coefficients a^"^ 7*^"^ and a*^"^ 

are not constrained by symmetries and must be determined from elsewhere. However for 
n = 1 they are fixed by the isospin charge of the hadrons to be 

a(i) = +1 , = +2 , = +1 , 7^'^ = -3 , ct(i) = . (58) 

Comparing the operators in the QCD and PQQCD chiral Lagrangians, one finds that the 
coefficients a^"'^ and both contribute to A^'^^ in the QCD limit of PQQCD, but away 
from this limit the operators are independent. 

At the one-loop level there are contributions from counterterms involving one insertion 
of the quark mass matrix mg. 



PQr){n),3 \ An) -^kji . 3 in R ., 



rdk 

+ (^3^)^ str (^+) 

b^s^ B^^^ Bijk str (r|+ v^^v^^...v^^ - traces , (59) 

where the coefficients bi^\...bg^^ are to be determined. The only constraint that exists on 
the 6-"^ is that the isospin charge of each nucleon is absolutely normalized, and thus b\^^ = 0. 
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The forward matrix elements of ^'^O^J^]',^ ,, in the z-th nucleon at the one-loop level can 
be written as 



/PQ(-^in),3 \ 



Pi 



(n) 



nl 



l-S 



Vi 



(n),0 



(n) , M,j I (n),l 

Pi uji + yj 7]l + yi 7]l 



{n),0 , (n),j . (n),l 

c] + Vj c\ + yi c\ 



Ui v^^v^^...v^^ Ui - traces , (60) 



where p-"^ are the LO contributions. The factor of 1 — 6^^ appears in the higher order 
corrections because the isovector charge is not renormalized. The diagrams shown in Fig. |^ 
give the leading non- analytic contributions to the wavefunction renormalization, Wi, the 
vertex contributions, rji""^'^ , that are independent of the charges of the ghost- and sea-quarks, 
and to the vertex contributions, ?7j-"'''"''', associated with the charges of the ghost- and valence- 
quarks. For the proton we find 



N 



M , % 



(a) 



N 




(c) 



M , % 



( 1 



N N 



N ^ N 

(b) 

^ M , Z 

N ^ N 



(d) 

( I 

V / N 



(e) 



(f) 



N 



M , 1 



N 



(g) 



N 



M , X 

\ N 



(h) 



FIG. 3. One-loop graphs that give contributions of the form ~ niQlogmQ to the matrix ele- 
ments of the isovector twist-2 operators in the nucleon. A solid, thick-solid and dashed line denote 
a 70-nucleon, AA-resonance, and a meson, respectively. The solid-squares denote an axial coupling 



given in eq.(3c), while the solid circle denotes an insertion of the strong two-pion-nucleon inter- 



action given in eg. (26 ). The crossed circle denotes an insertion of the tree-level matrix element of 
0^*^/12 ■••Mn- Diagrams (a) to (f) are vertex corrections, while diagrams (g) and (h) give rise to 
wavefunction renormalization. 



pin) ^ 1 ( 2a(«) - ) 
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+ gig A (2L 

+ ^ (-^un — SLurf + 2L;„ + 3Lid + 2Lju + ^Lj^ + 3-R^„,^„ + GB^^^rj^ + 3R 



'Vd,rid/ 



gi9A 



+4 



+2/3^ •* (-^'ud ~ -^ju ~ -^iu ~ Rvu,m ~ Rvu,vu) 
q;*'"^ (3L„y — GLucL — Ldd + 5Ljd + Lju + hLu + Liu 



+2/5^") (Ljo! - Lju + Lid - Liu - 2Ldd - Rvu,vu - '^Rvu,Vd ~ ^dm) 



.2 I - (n) . 



+ q9an I 7 



{Jdd Juu ^Jud 2Jjd 2Jid %)u,r]u ^m^Vd ^•^m'?. 



+ gQ;^"'* ( ^Luu — 4:Lud — Ldd — 5Lju — 5Liu + Ljd + Lid) 

+ 2/5*'"'' ( Luu + Lud — 2Ldd — Lju — Liu + 2Ljd + 2Lid) 



+ '1 

8 

1 



a 



(n) 



{Ljd + Lju — Lud — Luu) + /^'•"^ {'^Ljd + 2Lju — ^Lud — 2Luu) 



91n l^"-^ 



{Juu ~l~ 2Jud 2Jjd Jju) 



1 1 

+-q;^"^ ( - Lud + 5Lju + Ljd) + ( -Luu - 2Lud + Lju + 2Ljd) 

gW^ (Liu - Lud) + {Liu - L, 



{Lid + - Lud - Ldd) + {^Lid + 2Liu - ALdd - 2Lud) 



+ g^AJV ( 7 



(-'ud + 2Jdd — 2Jid — Jlu) 



+ ( -5Lud - Ldd + 5Liu + L^^) + ^/^("^ ( -2Ldd - Lud + Liu + 2Lid) 
6 

-rriu { 26S"^ + 56^"^ - + 2>ht^ + 36^"^ 



im, - ht^ + i^i") - 26i") + 265 

1 



(m^- - ruu) Vi 



in) 
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cr'" = {mi - ma) 



(n) 



(61) 



where, for the contributions from the 70 intermediate states, we have defined Lab = 
m^jlog (m^^//i^), and Rx,y — H{Lx, Ly, Lx)- For the contributions from the 44 interme- 
diate states, we have defined Jab — JiiTiab, ^, /j), with 



J(m, A, ii) = (m^ - 2A2) log ( ^ j + 2AV - log 



'A - VA2 -m^ + ie' 



A + VA2 - m2 + ie 



(62) 



and T^^j^ = T-C{Jx, Jy, Jx)- For the neutron we find 



.(n) 
rn 



'In 



- ( 2«(") - ) 



9i9a 



+ 



a^"^ {Lud — ^Ldd — Ljd — Lid — 4-Rj,„,r)rf — '^Rr]a,Vd) 



+ 



q;^'*-' {6Lud + Luu — 'iLdd — Ljd — 5Lju — Lid — 



lu 



^Linu,Vu ^Lir)u,'nd ^Liqd,Vd) 



+ 2/3^"^ i^Luu + Ljd — Lju + Lid ~ Liu + Ltr^u,Vu + '^Ltr^u,Vd "I" Ltr^d,Vd) 



'In 



M,l 

In 



( -5Ldrf + 4L„d + + 5Ljd + 5Lid - Lj„ - Li^) 
+-^(^'^'^^ ( — -^dd — Lud + 2-L„„ — 2Lju — 2Liu + Ljd + Lid) 



^4 



a 



(») 



(-f'jd + Lju — Lud ~ Luu) + 2/3''"^ (Ljd + 2Lj„ — Lud ~ ^Luu) 



1 



An) ' 



+7;9an \ T^""* ^ 1 {Jud + 2Juu - 2Jju — Jjd) 



+ la^^^ ( -5Lud - Luu + + L,u) + '-P^^^ ( - 2L„. + L,. + 2L,„) 



^1 



(L;<i + - Ldd - Lud) + 2/3^"^ {Lid + 2Liu - Ldd - 2Lud) 



(n) 
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1 / . ^ (t(")\ 

+^91n (7^"^ - ^ j {Jdd + '2Jud - 2Jiu - Jid) 
o o 



1 
3 

+i (m, + - 2h'^< 



26 



26(") + 36 



(n) 



(63) 



In the QCD and isospin limits, these expressions reduce down to those obtained in 
Refs. i25|J2r 



nl^ 



9 



27 



traces 



(64) 



where we have not shown the contribution from local counterterms involving a single in- 
sertion of rriq. Notice that the QCD limit is independent of the coupling g\, as required. 
Furthermore, in the isospin limit alone, where rnj = mi and m„ = rrid, these expressions 
reduce down to those obtained in three-flavor PQQCD p| when the strange quark mass is 
taken to be very heavy. 



VI. THE AXIAL- VECTOR CURRENT 



Matrix elements of the axial-vector current, 5"r"7^75g, are extensively studied on the 
lattice p4| . In xPT, there have been numerous computations of these matrix elements |ll^j32 



at the one-loop level, both including and excluding the decuplet as an intermediate state. 
In analogy to the extension of the electric- and isovector-charge matrices to PQQCD, there 
is a non-uniqueness in the extension of the isovector axial currents, Qt""^ fi'^^Q ■, to PQQCD. 
We use the charge matrix in eq. (|56| ) for the flavor- conserving currents, and for the flavor- 
changing current we replace with r+ in the upper 2x2 block. At LO in the chiral 
expansion, the axial current takes the form 



j;,5^ 2a (bS.B r|+) + 2(3 {bS, tI^B) + 2H {t'S, r\_,% 




>3 



+ .... 



(65) 



At one-loop level the matrix elements of the axial current between nucleons of flavor "a" 
and "b" are 
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Pab + 



167r2/2 



1 r -, (j) (I) 

Vab - Pab-[Wa + Wb] + Vj r]^h + Vl v j 



+ Cab + Vj cl^b + yi C^ab 



2UtSMa 



(66) 



where the wavefunction renormahzation contributions are given in eqs. (|6TD and (|63D. The 
Cab's correspond to contributions from local counterterms involving one insertion of the mass 
matrix, tuq, 



J fi,5 



S„B 



+ 5^ ^^y.(m+v,) b"'^ (^3^)^ (^^)« + 8"'' (r|+)' str {M^ 

+ 68 b"'' S^B,,k str (r|+ 



(67) 



where the coefficients, 6,-, must be determined from data or from lattice calculations. 



For the matrix element of between proton states, we find that the diagrams in Fig. 



give 



Ppp = Qa 



VuVd 



9A9l 
12 
9! 



10 



—9an9aa [ QJud + Juu - Jdd + 2Jjd + 2 Jjrf + 7"^^^^ + T^^ria - ^Xj^nd 



+ g9AN9A [ Kud + Kuu + Kju + -f^/u + 2Sr,^riu " 2iS^^ 



„(i) 

Ipp 



y 

S'A (-^ju ~ Luu) + ^ [ -^jd + Lju — Lud ^ Luu ] 



~'^9aa9an i^Jjd ~ 2Jud + Jju ~ Ju 

oi 



^9 A + ^1^1 + 



^1 

+ 1^ [ 5Luu + QLud — 5Lju — 9Ljd 
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(a) 



(b) 




(c) 



M , Z 




(d) 



M , X 



( 1 

N V ^ ✓ N 
■g ' 

(e) 



M , Z 



N 



(f) 



N 



N 



M , Z 
N 

\ N 



(g) 



FIG. 4. One-loop graphs that give contributions of the form ~ mg log mg to the matrix elements 
of the axial current in the nucleon. A solid, thick-solid and dashed line denote a 70-nucleon, 
44-resonance, and a meson, respectively. The solid-squares denote an axial coupling given in 
eq.^^S^, while the crossed circle denotes an insertion of the axial current operator. Diagrams (a) 
to (e) are vertex corrections, while diagrams (f) and (g) give rise to wavefunction renormalization. 



„(0 
Ipp 



^pp 



"pp 
pp 



2 



ud 



K 



9a {L 



lu 



L 



ud) 



9i 



+ — [ Lid + Liu 



2Kjd 

~ Lud ■ 



^dd 



— -^9AA9aN C^Jld — '^Jdd + Jlu — Jud) 
oi 



+ 2 {Lud — Liu) 



^9\ + 9i9\ + \9A9l 



9i 

+ ^ [ 5Lud + QLdd — 5Liu — 9Lid 



2 

■7,9an9i [ Kud + '^Kdd - Kiu - 2Kid 



'2 



b3 + h + hi 



/ 4, r 
+ rrid ( --61 - -62 



+ {rrij + mi) 

bs {nij - niu) 
bs {mi - md) 



(68) 



where we have defined the loop function, K, 



ab 



K{mab,A,fi) with 



20 



m 



-A^ lo, 



m 



-AVA^ — log 



'A - VA2 - nr- + ie' 
A + VA2 -m2 + ie. 



+ 



3X 



Trm — V A^ — log 



/ A- VA^-m^ + ie ^ 
VA + VA2-m2 + ie^ 



(69) 



and = HiKaa, Khh, Kx)- For the matrix element of between neutron states, we 

find 



Pnn — —9A 

Vnn — ~9a [ ^ud ~ L^d — R 



„(i) 

Inn 



VdVd 



9i9a 



2 [ 2I/dd — + Ljd + Lid + '^R'VuVd + "^^ndm 



9A9l 
12 



[ — 2Lud — Ljd + 3Lju — L^^ + 3Li„ + ?>Rq^nu + ^Rnum + ^^Vd^id 1 



+^ [ -^dd ~ 2L„(i — 3L„„ — Ljci + 3Lju — Lid + 3L;„ ] 
+ 9a [ Lud — Ldd + Ljd + Lid ] + ^ [ Luu — Ldd + Ljd — Lju + Lid — Liu 

+ -^9an9AA [ QJud + Jdd ~ Juu + 2 Jj„ + 2 + + %i^r)d ~ '^'^VuVd \ 



7:9an9a [ Kud + Kdd + Kjd + i^T^d + 2S, 
y 



~ g9AN9l [ ^^dd + Kud — '^Kuu + Kjd — 2Kju + Kid — 2Kiu + ^Sr^^Vd ~ 4»S^„^„ ] 

9a {Ljd — Lud) + ^ [ Ljd + Lju — Lud ~ Luu ] 



77« 
Inn 



w 

nn 



81 



9an9aa {Jjd Jud + 2 2 J„„) 



+ - {Lud - Ljd) 



\9a + 9i9a + 



+ fl [ ^Lud 



9Luu - 5Ljd - 9L 



~ q9an9i [ Kud + 2Kuu — Kjd — 2Kju ] 

9a {Lid — Ldd) + ^[ Lid + Liu — Ldd — Lud 



~'^9an9aa {Jid — Jdd + 2Jiu — 2 Jud) 
oi 

+ ^ {Ldd — Lid 



\9\ + 9i9\ + \9A9l 



9i 

+ [ 5Ldd + QLud — 5Lid — 9Liu 



2 

- q9an9i [ Kdd + 2Kud 



Kid - 2K 



lu 



4 112 2 



/2 5 1 
rrid (-W + - + h + h 



- {ruj + mi) ( -67 - -h 

bs {rrij - niu) 
bs {mi - md) 



(70) 
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For the np matrix element induced by the r"*" axial current, we find 

Pnp = 9 A 

Vnp = ~9A^riu'nd 

+ + ^ [ 2L„ci — 2Luu — 2Ldd — Lju — Ljd — Li^ — Li^ 

3 

+ ^ [ 2Lud + Luu + Ldd — Lju — Ljd — Liu — Lid ] ~ ^ t -^J" ~'~ ^^"^ ~^ ^i'^ ~^ "^'"^ ] 

^ '^9an9AA [ 4:Jud + Juu + Jdd + Jju + ^in + + Jld + 27^,^^^ + 27^^^^ — ATr/^r/d 
oi 

4 



+ n9AN9A [ 4:Kud + Kju + -ft'iu + -ft'jd + -ft'id + ^^^r?^ + ^r^dVd ~ ^'5, 



+ 7:9an9i [ 2-ft'ud + -R'mu + -ft'dd — -^jM — Klu — Kjd — Kid ] 
y 

(m„ + m^) (^-ifoi + ^62 - + + ^^5^ + {mj + nil) (^^br - ^bfi^ , (71) 

along with r^^^ = r^g = eg = eg = 0. 

These three independent matrix elements each reduce down to that of QCD when rrij — > 
TTiu and mi — > rrid. In the isospin limit and taking A — the axial-vector current matrix 



elements in QCD, at one-loop in the chiral expansion, are ||T4p2 

Ljj- f / 2 \ 2 2 ^0 2 A f \ 

Tnp = 9A - y9A\}^ 2g^ j + -gA9AN + -^9aa9an J , (72) 

where we have not shown the contribution from local counterterms involving one insertion 
of rUq. We have used the fact that K^, —>■ in the A ^ limit. 



VII. CONCLUSIONS 

The properties of the neutron and the proton have provided an important benchmark 
for comparing the predictions of continuum hadronic effective field theories (and models) 
with nature. They now provide a benchmark for the ever increasing number of lattice QCD 
calculations, be they quenched, partially-quenched or unquenched. In this paper we have 
analyzed several observables in two-fiavor PQ^PT to allow for extrapolations of current and 
future lattice simulations of two-fiavor QCD. One- loop level computations of the nucleon 
masses, magnetic moments, axial-vector currents, as well as the forward matrix elements of 
isovector twist-2 operators that are directly related to parton distribution functions, have 
been presented. 
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